Spatial distribution of persistent sites 



G.Manoj and P.Ray 

The Institute of Mathematical Sciences, C. I. T. Campus, Taramani, Madras 600 113, India 

(February 1, 2008) 



Abstract 



We study the distribution of persistent sites (sites unvisited by particles A) 
in one dimensional A + A — > reaction-diffusion model. We define the empty 
intervals as the separations between adjacent persistent sites, and study their 
size distribution n(k, t) as a function of interval length k and time t. The decay 
of persistence is the process of irreversible coalescence of these empty intervals, 
which we study analytically under the Independent Interval Approximation 
(IIA). Physical considerations suggest that the asymptotic solution is given 
by the dynamic scaling form n(k,t) = s~ 2 f(k/s) with the average interval 
size s ~ i 1 / 2 . We show under the IIA that the scaling function f(x) ~ x~ T as 
x — > and decays exponentially at large x. The exponent r is related to the 
persistence exponent through the scaling relation r = 2(1 — 8). We compare 
these predictions with the results of numerical simulations. We determine the 
two-point correlation function C(r,t) under the IIA. We find that for r <C s, 
C(r, t) ~ r~ a where a = 2 — r, in agreement with our earlier numerical results. 
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I. INTRODUCTION 



The persistence of fluctuations in stochastic processes has been an important topic of 
study in recent times Of primary interest in this context is the persistence probability 
P{t), which is the probability that a given stochastic variable <j)(t) with zero mean retains its 
sign during the time interval [0,t]. A power-law decay P(t) ~ t~ e is found in many systems 
of physical interest. Consequently, much effort has gone into the calculation of the new 
exponent 9 and studying its properties. There have also been several experimental studies 
of the persistence exponent 9 in real life systems 0. 

A particulary important class of systems whose persistence behaviour has been inves- 
tigated are spatially extended systems with a stochastic field 0(x, t) at each lattice site x. 
The time evolution of 0(x, t) is coupled to that of its neighbouring sites. 0(x, t) could be, 
for instance, an Ising spin @,||, a phase ordering field 0, a diffusing field or the height 
of a fluctuating interface J7||. The persistence probability P(t) is then the fraction of sites x 
where 0(x, t) has not flipped sign till time t. Recently it was observed that the set of per- 
sistent sites form a fractal and the time evolution of their spatial correlations obey dynamic 
scaling M. The purpose of this paper is the investigation of these spatial correlations. 

For concreteness, we study the one-dimensional A + A — > model. Our primary motiva- 
tion for this choice is the simplicity of the dynamics of the model which makes an analytic 
approach possible. In addition, this model is closely related to the d = 1 Glauber-Ising 
model, which is perhaps the only non-trivial model where 9 is known exactly. We study 
the distribution of the separations between nearest neigbour pairs of persistent sites. We 
call this the Empty Interval Distribution n(k,t), defined as the number of occurences where 
consecutive persistent sites are separated by distance k at time t. This quantity is a direct 
probe of spatial correlations in the distribution of particles. If the particles are distributed 
at random with some average density p, then n(k) = p 2 (l — p) k ~ e~ Afc at all k. Any slower 
mode of decay is indicative of spatial correlations. 

In this paper, we study the time evolution of the size distribution n(k, t) of these Empty 
Intervals. Persistence decay is identified with the irreversible coalescence of these intervals. 
The paper is organised as follows. In the next section we write a rate equation for the coa- 
lesence of these intervals under the approximation that the lengths of adjacent intervals are 
uncorrelated (HA). We give phenomenological arguments about the asymptotically relevant 
dynamical length scale as well as the coalescence probability. These arguments, combined 
with the rate equation gives the dynamic scaling behaviour of n(k,t) at late times t. We 
compare our predictions with numerical results. In section III, we use the IIA to predict the 
two-point correlations in the distribution of persistent sites. The predictions are found to 
be in agreement with recent numerical results, showing that IIA is valid. 

II. THE EMPTY INTERVAL DISTRIBUTION 

In the A + A — > model, a set of particles are distributed at random on the lattice with 
average density no- Over one time step, all the particles make an attempted jump to either 
of the neighbouring sites with some probability D. If two particles meet each other, both 
disappear from the lattice. In one dimension, the density of particles decay with time as 
n(t) ~ (87r.Dt)~2 as t —>■ oo [f|. Persistent sites in A + A — > model at any time t are defined 



2 



as the sites which remained unvisited by any diffusing particle throughout the time interval 
[0 : t}. Empty Intervals (which we will call 'Interval' for simplicity henceforth) are defined as 
the separations between two consecutive persistent sites. By definition, an Interval cannot 
contain a persistent site, although it may contain one or more diffusing particles A. The 
total number (per site) of Intervals of length k at time t is denoted by n(k, t) and is called 
the Empty Interval Distribution. 

To start with, the the particles are put randomly on the lattice so that n(k, t = 0) = 
tiq(1 — n ) k ~ e~ Xk where A = — log(l — n ). With time, the particles diffuse on the 
lattice, making the sites non-persistent. n(k, t) evolves satisfying the following normalisation 
conditions. If I m (t) = J2h k m n(k, t) ~ f£° n(s, t)s m ds is the m-th moment of the distribution, 
then 

I {t) = P(t) ~ r d ; h(t) = l I 2 (t) = s(t) (1) 

The first condition follows from the definition of n(k,t), the second one implies length 
conservation and the third condition gives the mean interval size s(t). The probability 
distribution of interval lengths is p(k,t) = ^ = P(t) -1 ri(A;, t) so that J2kP{k,t) = 1. 

Two neighbouring Intervals can coalesce when the persistent site between them is de- 
stroyed by a diffusing particle at the boundary of either of the Intervals. Note that this 
coalescence process is irreversible. For simplicity, we consider only binary coalescence in a 
single time step where two adjacent Intervals of lengths k% and k 2 , separated by a persistent 
site, coalesce and form a new Interval of length k\ + k 2 when the persistent site is 'killed' by 
a particle (Fig. I). To study this process analytically, we invoke a mean-field approximation 
- the lengths of adjacent Intervals are treated as uncorrelated random variables with proba- 
bility distribution p(k,t). This is the Independent Interval Approximation (HA), which has 
been used to study a variety of problems in one dimension [pT|,PH|. 



A. Rate Equation for Interval coalescence 

Assuming that IIA is valid, the time evolution of n(k, t) is given by the rate equation 

dn(k,t) 1 ^ , N „ w„ , 

- n{m, t)p{k — m, tjK [m, k — m,t) — 



dt 2 m=1 



n 



(k,t)Y,P(m^)K(m,k,t) (2) 



m=l 



where i^(m 1 ,m 2 ,t) is the probability that two adjacent Intervals of lengths mi and m 2 
coalesce at time t. The first term in Eq.0 represents the increase in number of Intervals 
of size k through coalescence of smaller intervals, while the second term is the loss term 
representing the decrease in number when Intervals of size k merge with other Intervals. 

To solve the above equation for n(k,t), one need to know the form of the reaction 
kernel A'(mi,m 2 ,t). The process of coalescence of Intervals involves the destruction of the 
persistent site in between them by a particle, which can come from either of the Intervals. 
So, quite generally, 

K(m 1 ,m 2 ,t) = Q(m 1 ,t) + Q(m 2 ,t) (3) 
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where Q(m,t) is the fraction of intervals of size m which is destroyed at time t. Q(m,t) 
satisfies the following condition by definition. 



^2n(m,i)Q(m,t) 



dP(t) e 



dt 



P(t) 



(4) 



where we have made use of the fact that P(t) ~ t~ e . 

The form of Q(m, t) can be argued for in the following way. An Interval of length m at 
time t can contain a particle anywhere inside it only if the interval length is at least of the 
order of the diffusive scale \fDt. That is, Q(m,t) ~ for m <C \f~Dt. It is also known that 
the particle distribution is correlated over length scales r <C wheras it is completely 

random over r \f~Dt |T|. So we expect that for m ^> \^Dt, Q(m, t) — > ct(t), independent 
of m. These physical considerations leads us to suggest the following dynamic scaling form 
for Q(m, t). 



Q(m,t) = a(t)/3( 



rn 



Dt 



where the function (3(x) is expected to have a sigmoidal form, ie., /3(x) 
and j3{x) — > 1 for x ^> 1. The function a(t) will be determined later. 



(5) 



for x -C 1 



B. Dynamic scaling 

We assume that at asymptotic times, the distribution n(k, t) is characterised by a single 
dynamic length scale s(t). We note that there are two relevant length scales in the problem. 
The first is the diffusive scale £_d(^) ~ \f~D~i entering the scaling form Eq.|5|for the coalescence 
probability. On the other hand, the inverse of the persistent fraction P(t) is also a length 
scale, which we shall call the persistence scale, denoted by C p {t) ~ t e . The asymptotic 
behaviour is expected to be dominated by the larger of the two, ie., the diffusive scale Cjy(t) 
in the present case (since 9 < 1/2). 

We now invoke the dynamic scaling ansatz, ie., n(k,t) oc /(-) with 



.l/z 



(6) 



From the length conservation condition given by the second part of Eq. [T] it follows that 
the prefactor is ~ s~ 2 . Thus, the scaling solution for n(k,t) is written in the form 

n(k,t)=s(ty*f(J^j (7) 

Substituting Eq. |7| and Eq. |5| in Eq. [|, we find a(t). 

»W - ^ (8) 

where B = J °° (3(x)f{x)dx. Substituting Eq.[7| in the normalisation conditions Eq.[T|, we 
find the following conditions on the scaling function. 
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f(x)dx = sP(t) 



f(x)xdx 



(9) 



In the first integral, the lower limit is set as s(t) _1 to take care of any possible small 
argument divergence. 

Substituting Eq.|| [7|, || and || in Eq. |, we find the following equation for the scaling 
function f(x). 



■qdf_ 
z drj 



e 

B J s (t) 



f{%)f{v ~ x)[j3{x) + P(ti - x)]dx 



--0- ^s{t)P{t)(5{rj) 
z B 



(10) 



where the scaling variable 77 = ^£y. 



Case I: r\ >C 1. 

For r] <C 1, all /3(x) ~ for x < rj. This case corresponds to small Intervals, ie., those 
which are not large enough to contain a diffusing particle till time t. In this case, the 
equation reduces to 77S = —(2 — z8)f(r]) which has the solution f{rf) ~ v^ T where the new 
exponent r is related to 9 through the scaling relation 



r = 2 - zO 



(11) 



For the model under consid- 
5/4. 



From Eq. |7| this implies that for k s, n(k,t) ~ t 9 k 
eration here, 9 is known exactly to be 3/8 0] which gives r 
Case II: r] > 1. 

For general values of r), f3(rj) is non-zero, and because r > 1, the first integral diverges 
near x = as x~^ r_1 ' ) . There is another divergence in the last term, of the form t 1 ^^ 9 . 
It can be shown that this term can be exactly cancelled by the divergent part of the first 
integral. After carrying out this 'regularisation' (details to be found in Appendix A) and 
putting z = 2 in Eq.[10] the equation for the scaling function f(r)) stands as 



2 T = ^B Jo f{x)fiv ~ X)[m + ^ ~ X) ~ m]dX 



^(3( V ) f*f{x)\f{ri-x)-f(rj)]dx 
B Jo 
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B' 



P(v) / f{x)dx 



(12) 



A general solution of this equation requires the knowledge of the detailed form of the 
scaling function f3{r]). However, for large values of rj where (5{rj) ~ 1, one can simplify this 
equation. We define the point 77* sufficiently large such that for 77 > 77*, /3(x) = 1 within the 
limits of accuracy required. Without any loss of generality, one can put 77* = 1 by rescaling 
the length scale s(t) accordingly. For 77 > 1, we define f(r)) = h(r]), whose equation is 



rj dh 
29t7 



9 
B 



2 I" h(x)h(r) 



x)dx + / f(x)[h(rj 
Jo 



x 



h(rj)]dx 



[l-20)h{v) (13) 



This equation has a solution of the form h(i]) = Ge Xv as can be shown by direct 
substitution. The constants G and A are related through the relations 
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XB = 26G 



(14) 



and 



\ + 26 = 1 + -F(A) 
B 



(15) 



where F(X) = $ f(x) \e Xx (l + (3{x)) - 1 



dx and 



B 



1 G 

f(x)p(x)dx + — e" 
o A 



(16) 



by definition. Eq. [T4] - |lT)| formally gives the constants A and G. However, the actual 
evaluation of these constants requires the knowledge of the function f(x) in the entire range 
[0:1] (and not just near x = 0, where f(x) ~ x~ T ), which, in turn, is possible only if 
the detailed form of f3(x) is known. Hence we will restrict ourselves to showing that the 
parameter A > 0, which is required for the solution to be physically reasonable. 

In Eq. [H], we note that B > ye _A , depending on how sharply f3(x) rises near x = 1. 
The equality holds for the step function f3(x) = Q(x — 1) where Q(x) = for x < and 
Q(x) = 1 for x > 1. After using this inequality in Eq.|14|, we find that A > — log(2#). Since 
9 < 1/2, it follows that A > 0. 



C. Numerical Results 

We determine the distribution n(k, t) numerically by simulating v4 + A — > model on one 
dimensional lattice of size N = 10 5 with periodic boundary condition. Particles are initially 
distributed at random on the lattice with some average density no, and their positions 
are sequentially updated — each particle is made to move one step in either direction with 
probability D = 1/2. When two particles meet each other, both are removed from the 
lattice. The time evolution is observed up to 10 5 Monte-Carlo steps (1 MC step is counted 
after all the particles in the lattice were touched once). The simulation is repeated for 
several random starting configurations of the particles for any particular initial density and 
we repeat the entire simulation for four different initial density uq. For any uq, we determine 
the number of intervals of length k (per site) at time t. 

To compute the mean interval size s(t), we ran the simulation upto t = 10 5 time steps, 
and averaged the results over 100 starting distributions of particles, with the same initial 
density. In Fig. II, we plot s(t) vs t for four different values of n - 0.2, 0.5, 0.8 and 0.95. For 
hq = 1/2, we find that s(t) ~ at 1 ^ with z ~ 1.97(1) and a ~ 5.96. But for other values of 
no, we find that the observed value of z is different from 2. In Fig. Ill, the running exponent 
d(logs) /d(logt) is plotted with l/(logt) and the results show the systematic deviation away 
from the value 1/2 expected from the scaling picture presented in the previous section. We 
will discuss about the possible origin of this deviation later. 

In Fig. IV, we plot the scaling function f(x) = s(t) 2 n(k,t) against the scaling variable 
x = k/s(t) for t = 10 4 and 3.10 4 . To get the nature of the scaling function one needs to 
average over lots of configurations. This has restricted us to smaller time steps and data 
for three values of initial density, n = 0.2, 0.5 and 0.8 averaged over 500, 1000 and 1500 
different initial distribution of particles respectively. For all no, we find that the scaling 
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function f(x) ~ x~ T for x « 1 and decays exponentially for higher values of x. For 
n = 0.5, we find r = 1.25(1) in accordance with the scaling relation Eq |ll|. For uq = 0.2, 
we find r ~ 1.32(2) while for n = 0.8, the observed value of r is 1.13(2). For all n , r 
satisfies the scaling relation Eq. [TT] if z is replaced by its effective value. 

For general values of no, we find that the numerical values of s(t) supports the following 
form (within the time range studied) 

Sno (t)~a^ + 6K)^. (17) 

The non-universal constant b is <,= or > for <, = or > 0.5. To compute the 
prefactor b and the exponent (ft, we plot the difference As no (t) = \s no (t) — si/ 2 (t)| vs t, for 
n = 0.2,0.8 and 0.95 (Fig. V). The exponent cf) is numerically found to be close to the 
persistence exponent 6 = 0.375 (Table I). We find that as n — > 1, the constant b undergoes 
a sharp rise so that the effective dynamical exponent of s(t) is numerically close to 9 for an 
appreciable range in time (Fig. III). At the same time, we note that only the first term in 
Eq.|T^ is asymptotically relevant since <t> < \- 

The two terms in Eq.pj] can have their origin from the two dynamical length scales in 
the problem, the diffusive scale £_d(£) ~ t 1 ^ 2 and the persistence scale C p (t) ~ t e . For large 
no, the typical interval length between two consecutive persistent sites is determined by the 
decay of persistence only, rather than the diffusion of the particles. So, it is understand- 
able that the dynamical behavior of s(t) coincides with that of C p (t) at least at the initial 
times. However at late times, when the particle density falls down as a result of annihila- 
tion, the situation becomes same as that of starting with low no and the decisive scale is 
£i)(t). However, the precise form and behavior of the prefactor b(no) with no remains to be 
understood. 



III. TWO-POINT CORRELATIONS 

A good picture of the spatial distribution of the persistent sites and the presence of any 
possible correlation in their distribution is obtained from the two-point correlation C(r,t), 
which is defined as the probability that site x+r is persistent, given that site x is persistent 
(averaged over x). 

C(r,t) = (p(x,t))- 1 (p(x,t)p(x + r,t)) (18) 

where the brackets denote average over the entire lattice and p(x, t) is the density of 
persistent sites: ie., p(x, t) = 1 if site x is persistent at time t, and otherwise. Clearly, 
(p(x,t)) = P(t) by definition. 

Within the IIA, the relation between C(r,t) and n(r,t) (We consider r ^> 1, so that the 
discreteness of the underlying lattice can be ignored) can be written as the following infinite 
series: 

C(r, t) = P(t) _1 n(r, t) + P(t)~ 2 £ dx n(x, t)n(r -x,t) + 

/r pr—x 
dxn(x,t)J dyn(y,t)n(r-x-y,t) + .... (19) 
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The first term corresponds to the case where there is no other persistent site in the range 
[0 : r], ie., a single Interval of length r. The second term gives the probability that the range 
is split into two Intervals of length x and r — x by the presence of a persistent site at x, the 
third term gives the probability for three Intervals and so on. 

The above series can be rewritten as the following self-consistent equation for C(r,t). 

P{t)C{r, t) = n(r, t) + J n(x, t)C{r - x, t)dx (20) 

In terms of the Laplace transforms C(p,t) = C(r,t)e~ pr dr and n{p,t) = 
n(s,t)e~ ps ds Eq. p0| becomes 

From Eq.fj], we find 

n(p,t) = s- 1 f(ps) (22) 

where f(q) = J s °°i f(r])e~ qT] dr). which can be written in the following regularised form, 
using Eq.^|. 

f(q) = s(t)P(t)-f 1 (q) (23) 

where 

/•oo 

fM= m[l-e- qr >]d V . (24) 
jo 



Substituting Eq.E2l E3 and Eq.E3 into Eq.ET| we find that 



The second term in RHS can be neglected at late times, since s(t)P(t) diverges as t 1 ^^ 9 . 
It follows that in this limit, C(r,t) has the dynamic scaling form 

C(r,t) = P{t)g U-yj (26) 

where 

m ^ m (27) 

is the Laplace transform of g(x): g(q) = / °° g(x)e~ qx dx. 

The preceding expressions can be used to deduce the limiting behaviour of the scaling 
function g(rj) for the cases rj 1 and rj ^> 1, without needing to solve Eq.[L9| or ^0| explicitly. 
Case I: rj 3> 1. 
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To find the asymptotic behaviour of g(r]), we note that fi(q) vanishes near q = as 
fi{q) ~ <?• Thus g(q) ~ - as q — > from Eq.[27]. By standard results in the theory of Laplace 



transforms this implies that g(rj) ~ 1 as 77 — ► 00. 
Case II: 7? < 1. 

To analyse this case, consider the real-space relation Eq.[19|. For 77 <C 1, or equivalently, 
r < s, we have shown that n(r,t) ~ P(t)r~ T . It is clear that in this range, the RHS of 
Eq.[19] is time independent, so C(r,t) in the LHS should also be time independent. From 
the dynamic scaling form Eq.|26], we find that this is possible only if the scaling function is a 
power-law near the origin: g{rj) ~ r/ _a as 77 — > 0. After substituting in Eq.^ and requiring 
the resulting expression to be time independent, we find 

a = z6 (28) 

We find C(r,t) ~ r~ a for rCs and C(r,t) ~ P(t) for r ^> s. The power law decay at 
small distances is expected, because the RHS of Eq.|l^ contains only scale invariant terms 
in this limit, hence the LHS also should be likewise. In Appendix B, we show that this is 
also consistent with Eq ,p0|. 

We see that in the IIA calculation, the length scale s(t) demarcates the correlated and 
uncorrelated regions for C(r, t). In the correlated region (r <C s(t)), the persistent sites form 
a fractal with fractal dimension df = d — a = 7, with the correlation length s(t) increasing 
with time as s ~ t 1/2 . The IIA results agree very well with that of numerical simulations 
||, showing the validity of the approximation. 



IV. CONCLUSION 



Persistent sites are shown to have strong correlations in their spatial distribution. In 
one dimensional A + A — > reaction-diffusion system, we show that there is a length scale 
s(t), diverging with time as s(t) ~ t 1 ^, which demarcates the correlated region from the 
uncorrelated one. We argue that z = 2 at large t limit. Persistent sites separated by 
distance < s(t) are highly unlikely to have a particle A between them and so retains their 
persistent character. Only persistent sites separated by distance ^> s(t) take part in the 
decay of persistence at subsequent times. 

We find that if k is the distance of separation between any two consecutive persistent 
sites, then for k s(t), the distribution of k is scale-free and decays algebraically as k~ T 
with t = 2 — z9. We show this using the IIA (Independent Interval Approximation), which 
assumes no correlation in the lengths of any two adjacent intervals. We have verified our 
results by numerical simulations which suggests the validity of the IIA. Under the IIA, our 
calculation for the two-point correlation shows that over length scales r s(t), the persistent 
site distribution over the lattice is a fractal with dimension df = r — 1, in accordance to our 
earlier observations ||. 
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APPENDIX A: 



The divergence in the first integral in Eq. |TI] can be separated out as follows. We 
write f(ij - x) = f(rj) + A x f(rj) and P(rj - x) = P(rf) + A x f3{r]). so that lim x ^ A x f(rj) = 
lim x ^ A x p(r]) = 0. 

After substituting for f(r) — x) and (3(t] — x), the divergent part of the integral separates 
into the following terms. 

I" f(x)f{v - x)(3(v ~ x)dx = f( V )P{v) I" f(x)dx + f( V ) I" f{x)A x (3{r,)dx + 

Jsit)' 1 Js(t)- 1 JO 

(3( V ) [* f(x)A x f( v )dx+ \ l f(x)A x f(rj)A x p( V )dx 
Jo Jo 

The first term is divergent near the origin, while all other terms are finite by construction. 

Now we rewrite the first term using the equality /jT^-i f(x)dx = s(t)P(t). After some 

simplifications, the integral becomes 

P f(x)f(v ~ x)/3tn - x)dx = f(r))(3(r))s(t)P(t) + p f(x)f( V - x) \J3{ V - x) - (3( V )\ dx + 

Jsit)- 1 Jo 

rj 

JJ f(x) [f(v -x)- f(v)\ - (3(v)f(v) J™ f(x)dx 

The first term is the divergent part of the integral, which exactly cancels the last term 
in Eq. [10|, to give the regularised Eq. |H| 



APPENDIX B: 



For r 3> 1, it is reasonable to assume that the higher order terms in the RHS of Eq.[l9| 
will contribute more than the first term, ie., the range [0 : r] is more likely to be covered 
with more than one Interval than a single one of length r. After using this approximation, 
and substituting n{r,t) ~ (r — \)P{t)r~ T in the continuum limit, Eq. ^ is simplified to 



C(r,t)~(r-1) f V -x)- T C(x,t)dx 



Our purpose is to see if the equation 



r-l 



(r- 1) / x~ a {r-x)- T dx (Bl) 



is consistent for a = z9 = 2 — r (Eq.|TTD at r ^> 1. 

The integral I = f[~ x~ T (r — x)~ a dx can be transformed by change of variables into the 
more standard form 



rr-2 r l-a 

/ (1 + yY T [r - 1 - yY a dy ~ F(1,t; 2 - a; — r) for a < 1 and r > 1. 

Jo 1 — a 

where F(a, b; c; z) is the Gauss Hypergeometric function. For b = c, F(a, b; b; z) 
(1 — z)" a exactly, independent of b |14[]. Thus, for a = 2 — r we find 

(r - 1)1 = r-«[l + o(~)] 
which is consistent with Eq. plj , at r 3> 1. 
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TABLES 



no 


b 




0.20 


-6.621 


0.34372(11) 


0.80 


15.701 


0.35495(5) 


0.95 


84.672 


0.36572(4) 



TABLE I. Results for the prefactor b and exponent (j) as measured from simulations. The 
numerical value of 4> is found to be close to the persistence exponent 9 whose exact value is 0.375. 
The figures in brackets represent statistical error in the last decimal place. Note the sharp rise in 
b as no — > 1. 
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FIGURES 
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FIG. 1. In the picture, white circles are persistent sites (numbered 1, 2, 3..) and dark triangles 
are diffusing particles. Two Empty Intervals £12 and E23 are shown to merge together to give a 
new Interval E13 when the persistent site 2 at the boundary is killed by a diffusing particle. 
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FIG. 2. The length scale s(t) is plotted as a function of time t. The straight line is a fit, with 
slope 1/2. 
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FIG. 3. The effective exponent d °^ t is plotte against 1/logi for four values of starting density. 
For no = 0.5, the exponent value is close to 0.5, expected from the scaling arguments. For other 
values of no, systematic deviations away from 0.5 is observed. 




FIG. 4. The scaling function f(i]) = s(t) 2 n(k,t) is plotted against the scaling variable 
r] = k/s(t) on a logarithmic scale. There is a power-law divergence at small rj and exponen- 
tial decay at large rj, as predicted by the IIA calculation. The observed value of r for uq = 0.8 is 
seen to be appreciably different from that for other no- 
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FIG. 5. The difference As(t) = \s(t) - s 1/2 (t)\ is plotted against t for n = 0.2,0.8 and 0.95. 
The straight line is a fit with slope 3/8. 
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